Superfluidity of spin-polarized 6 Li 
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We study the prospects for observing superfluidity in a spin-polarized atomic gas of 6 Li atoms, 
using state-of-the-art interatomic potentials. We determine the spinodal line and show that a BCS 
transition to the superfluid state can indeed occur in the (meta)stable region of the phase diagram if 
the densities are sufficiently low. Moreover, for a total density of fO 12 cm -3 , which still fulfills this 
requirement, we find a critical temperature of only 29 nK. We also discuss the stability of the gas 
due to exchange and dipolar relaxation and conclude that the prospects for observing superfluidity 
in a magnetically trapped atomic 6 Li gas are particularly promising for magnetic bias fields larger 
than 10 T. 
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Ultracold atomic gases have received much attention 
in recent years, because of their novel properties. For 
instance, these gases are well suited for high-precision 
measurements of single-atom properties and for the ob- 
servation of collisional and optical phenomena that reflect 
the (Bose or Fermi) statistics of the constituent particles. 
Moreover, a large variety of experimental techniques are 
available to manipulate the atomic gas samples by means 
of electromagnetic fields |lj , which offers the exciting pos- 
sibility to achieve the required conditions for quantum 
degeneracy and to study macroscopic quantum effects in 
their purest form. 

At present, most experimental attempts towards quan- 
tum degeneracy have been performed with bosonic gases 
and have been aimed at the achievement of Bose-Einstein 
condensation. In particular, most of the earlier experi- 
ments used atomic hydrogen These experiments 
provided crucial ingredients for the recent attempts with 
alkali vapors, for which the experimental advances to- 
wards the degeneracy regime were so rapid that Bose- 
Einstein condensation has actually been reported now 
for the isotopes 87 Rb @ and 7 Li g. 

In view of these exciting developments it seems timely 
to investigate theoretically also the properties of spin- 
polarized atomic 6 Li, since 6 Li is a stable fermionic iso- 
tope of lithium that can be trapped and cooled in much 
the same way as its bosonic counterpart. Therefore, mag- 
netically trapped 6 Li promises to be an ideal system to 
study degeneracy effects in a weakly-interacting Fermi 
gas, thus providing valuable complementary information 
on the workings of quantum mechanics at the macro- 
scopic level. Moreover, using a combination of theoret- 
ical and experimental Q results, accurate knowl- 
edge of the interparticle (singlet and triplet) potential 
curves of lithium have recently been obtained which lead 
to the prediction of a large and negative s-wave scatter- 
ing length a of —4.6 • 10 3 ao (do is the Bohr radius) for a 
spin-polarized 6 Li gas. 



This is important for two reasons: First, the fact 
that the scattering length is negative implies that at 
the low temperatures of interest (A 3> ry, where A = 
{2itTi 2 /raksT) 1 / 2 is the thermal de Broglie wavelength 
of the atoms and ry is the range of the interaction) the 
effective interaction between the lithium atoms is attrac- 
tive and we expect a BCS-like phase transition to a su- 
perfluid state at a critical temperature 
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with ep — h 2 k F /2m the Fermi energy of the gas [p|. 
Secondly, we see that the critical temperature depends 
exponentially on 1/A;f|a| which usually, when the mag- 
nitude of a is of the order of the range of the interaction 
ry, is very large in the dilute limit kpry <C 1. There- 
fore, it was previously concluded that the BCS transition 
in a dilute fermionic system (in particular spin-polarized 
deuterium) is experimentally unattainable |10|. However, 
with the anomalously large scattering length of 6 Li this 
conclusion needs revision as we will see now in more de- 
tail. 

Since 6 Li has an electron spin of s = 1/2 and a nuclear 
spin of i = 1, the Is groundstate of 6 Li consists of six 
hyperfine levels which are labeled by |1) through |6) in 
such a way that their energy increases at small magnetic 
fields. At zero magnetic field these correspond exactly 
to the states |/,m/) with a total spin / of either 1/2 or 
3/2 and a hyperfine splitting of 3a^//2 ~ 10.95 mK. In 
a conventional trapping experiment one tries to achieve 
both electron and nuclear spin-polarization by trapping 
only atoms in the doubly spin-polarized state |6) = \m s — 
1/2, m, = 1). In the case of 6 Li, however, such a proce- 
dure creates a gas in which the atoms only interact ex- 
tremely weakly because the Pauli principle now forbids 
s-wave scattering. Therefore, we cannot take advantage 
of the large negative value of the triplet scattering length 
and the BCS-transition temperature is unattainable. 
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To avoid this problem we need to trap two hyperfine 
states. This can be achieved most easily by using rela- 
tively large magnetic fields B ^> a^j / 'fj, e ~ 0.011T (fx e 
is the electron magnetron), because then the electron 
and nuclear spins are almost decoupled and the states 
|4) ~ \m, = 1/2, mi = -1), |5) ~ \m s = 1/2, rm = 0), 
and 1 6) can all be trapped. Moreover, the trap can, for ex- 
ample, be loaded by creating first a doubly spin-polarized 
gas in the conventional manner and subsequently apply- 
ing a microwave pulse to populate one of the other trap- 
ping states. Notice that because the gas is now only 
electron (and not nuclear) spin polarized, there is a large 
cross-section Aira 2 for the thermalizing collisions required 
for evaporative cooling. In addition, this implies that 
the gas is in thermal equilibrium in the spatial degrees 
of freedom, even though the spin degrees of freedom are 
not. Notice also that the electron spin polarization is 
not complete because the states |4) and |5) have a small 
admixture (of 0(ahf / fj, e B)) of |m s = — 1/2, = 0) and 
\m s = — 1/2, mj = 1), respectively. Therefore, two atoms 
do not interact solely via the triplet interaction. When 
we consider the lifetime of the gas, however, we find that 
the influence of this can be neglected if the magnetic field 
is larger than 1 T. 

Although it is possible to study any combination of the 
three trapping states, we will consider here only a gas in 
which the atoms are in a mixture of state |5) or state |6) 
because this minimizes the number of decay processes. 
Furthermore, we analyze here first the homogeneous case. 
The influence of the trapping potential will be discussed 
in a separate publication. Taking only s-wave scattering 
into account and following Ref. (ll]] to include all two- 
body processes, we can determine the thermodynamic 
properties of this gas by considering the hamiltonian 
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Here the field ipa(x) (ip a (x)) creates (annihilates) an 
atom at position x in a spin state |a) that has a (renor- 
malized) chemical potential [i' a and an average density 
n a . In addition, Vb = J dx Vt (x) is the zero mo- 
mentum component of the triplet potential and An = 
Vo(ip5(x)ip6(x)) is the equilibrium value of the appropri- 
ate BCS order parameter. 

A first condition for the observability of the BCS tran- 
sition is that it takes place in the (meta)stable region of 
the phase diagram pi] ]. The requirements on the densi- 
ties in the two hyperfine levels that result from this condi- 
tion of mechanical stability of the gas can be obtained as 
follows. For the normal state (where An = 0) the hamil- 
tonian is diagonal and we can immediately compute the 



pressure p of the gas by evaluating the thermodynamic 
potential density ft/V — —p. This results in 
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where e(k) = h 2 k 2 /2m is the kinetic energy of the atoms 
and P = l/ksT. Furthermore, the density of lithium 
atoms in the two hyperfine states is determined by 
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N a (k) , 
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introducing the notation N a (k) for the occupation num- 
bers which in our case are equal to the Fermi distribution 
function [e l3x + 1) _1 evaluated at e(k) — fj! a . 

In the degenerate regime, neglecting corrections of 
0((fcsT '/ 'ep) 2 ), the equation of state Eq. (0) can be easily 
inverted and we find for the pressure 
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For the mechanical stability of the gas we must require 
that dp/dn a > 0. Using the above result this leads to 
two conditions on the densities in the two hyperfine lev- 
els, namely 

715 - 6R ( 3n6 /l) and n& ~ 6R ( 3n5 /f ) ' 

which have to be fulfilled simultaneously. The line where 
one of the equalities holds is called the spinodal line and 
it is shown for 6 Li in Fig. |l|. Note that for the highest 
metastable densities we have fc^M = tt/2. Therefore, 
the ratio ksTc/eF is at most 0.23 and our determination 
of the spinodal line is self consistent for the temperatures 
of interest. 

Within this density region (i.e. also for ^ tiq) we 
can now consider the critical temperature of the gas. For 
that we need to derive the BCS gap equation. This can 
be achieved most easily by diagonalizing the hamiltonian 
by means of a Bogoliubov transformation and then cal- 
culating the equilibrium value of An = Vo(ip5(x)tpe(x)} ■ 
In the limit of vanishing An, or equivalently T —> T c , this 
procedure leads to the linearized BCS gap equation 
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dk 1 - N 5 (k) - N 6 (k) 
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where we have introduced the appropriate Fermi en- 
ergy €f — (Ms + /Lt@)/2. As expected, the integral in 
the left-hand side of Eq. (^) has an ultraviolet diver- 
gence due to the neglect of the momentum dependence of 
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the triplet potential in the hamiltonian. However, from 
the Lippmann-Schwinger equation for the two-body T- 
matrix [Q we find that this divergence is cancelled by 
a renormalization of l/V to 1/T 2B (0, 0; 0) = m/Airaft 2 . 
Therefore, the critical temperature is determined by the 
condition 
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which is free of divergences because of the Fermi distibu- 
tions in the Cauchy principle-value integral. 

We have solved Eqs. (0) and (^) numerically for a given 
density in the two hyperfine states. Furthermore, we 
have verified that the energy dependence of the triplet 
s-wave cross section is unimportant even though the scat- 
tering length is large and negative. Physically the latter 
is a result of the presence of an almost bound state near 
the continuum treshold and one might therefore have ex- 
pected a relatively strong energy dependence. Our final 
results are summarized in Fig. |l] where we have plotted 
the critical values of n 5 and Uq for three different temper- 
atures. (They can, of course, be calculated directly from 
Eq. (|l|) if tt-5 = fig.) In particular, we find that non-zero 
critical temperatures are only possible if the 'polariza- 
tion' \n$ — n§\/(nQ + n^) is less than 3fcsT c /2ei?. We 
also find for Scf = |a4 — Msl around 2ksT c a reentrance 
behaviour, because at a fixed value of £f the function 
T c (5£f) is multivalued. This is shown more explicitly 
in the inset of Fig. [l] and can be proven analytically by 
showing that the linearized gap equation reduces to 
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for ksT c <C Sep <C ef and thus explaining the square- 
root behavior of T c (<5£f) found in this regime. 

Finally, we also need to consider the lifetime of spin- 
polarized 6 Li. It is well-known that spin-polarized atomic 
gases are not completely stable because collisions be- 
tween the atoms can induce transitions from one hyper- 
fine level \a) to another. In the case of interest these 
transitions can be caused both by the central (singlet and 
triplet) interactions V c {r) and by the magnetic dipole- 
dipole interactions V d (r). In general, the zero tempera- 
ture rate constant of a transition a/3 — > a' (3' equals 
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with Pa>[]' the momentum, and t' and ml the angular 
momentum quantum numbers in the final state p3[ . 

The central interaction V c (r) cannot flip the electron 
or nuclear spin. Therefore, the only transition due to 
this interaction that is possible in a gas consisting of 
atoms in the hyperfine states |5) and |6) is 65 — > 61. The 
coupling matrix element between these states is propor- 
tional to the exchange potential V ex (r) = Vr{r) — Vs(r), 



i.e. the difference between the triplet and singlet po- 
tentials. More precisely, if B 3> a^f/fJ-e it is equal 
to -(a hf /n e B)V ex (r)/2V2 . Therefore, at sufficiently 
high magnetic fields the coupling becomes small and the 
T-matrix element in Eq. (|9|) can be calculated in the 
distorted-wave Born approximation leading to 
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where l^^s) are the exact scattering states of two 
lithium atoms with angular momentum quantum num- 
bers i and to, total electron spin S, and at the energies 
of the initial and final state, respectively. 




FIG. 1. The required conditions for the achievement of the 
BCS transition in the (ns, ne) plane. The dashed curve shows 
the spinodal line, and curves 1, 2 and 3 give the conditions at 
which the critical temperature is 0, 29 and 400 nK, respec- 
tively. The inset shows the critical temperature as a function 
of SeF for ep = 385 nK. 

The decay rate as a function of magnetic field follow- 
ing from this expression is given in Fig. |^. We give the 
rate constant only for magnetic fields larger than 1 T be- 
cause we have verified numerically that only then can the 
exchange interaction be treated in first order perturba- 
tion theory. Note that the rate constant is much larger 
than the rate constant for the similar process in magnet- 
ically trapped atomic hydrogen. This is again a result 
of the fact that in the initial (triplet) state there is an 
almost bound state near the continuum threshold which 
enormously enhances the amplitude of the triplet wave 
function. Clearly, this is the drawback of having such a 
large scattering length. Nevertheless, we see that a mag- 
netic field of about 10 T suppresses the decay of the gas 
to the extent that for a density of 10 12 cto~ 3 , the lifetime 
is of the order of seconds. 

Of the various magnetic dipole-dipole interactions the 
electron-electron interaction is most important. Never- 
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theless, it is so weak that it can always be treated by 
first-order perturbation theory Thus, if B ^> auf/^e 
we can neglect the hyperfine coupling between the elec- 
tron and nuclear spins and the total dipolar rate consist 
of the sum of a one spin-flip and a two spin-flip con- 
tribution, i.e. G d = G lsf + G 2sf . Due to the different 
spin-matrix elements and the difference in the energy re- 
leased in the transition, we find the convenient relation 
G 2sf {B) = 2G lsf (2B). Hence, we need to consider only 
the one spin- flip process, which we again treat in the 
distorted-wave Born approximation. For magnetic fields 
larger than 1 T this leads to 

nisi 6 R u m (^oMe) 2 I ,2 

G J ~—^2mii e B — A |r 20 | , (11) 

07T ft 

where r 2 o — J °° dr ip^i ' i r ) r ~ i i } oi > ( r ) 1S the relevant 
transition matrix element and the radial wavefunctions 
ip^\r) are normalized as 

- J^^JY^t) . (12) 
V irn r 

The rate constant for the one spin-flip process is also 
shown in Fig. || Again, it is much larger than the rate 
constant of a similar proces in atomic hydrogen. Notice 
that for magnetic fields larger than 10 T the electron- 
electron dipolar rate dominates the decay and then also 
leads to a lifetime of the order of seconds for a den- 
sity of 10 12 cm~ 3 . The decay rates due to the electron- 
nucleus dipolar interaction are even smaller by a factor of 
(/xjv//x e ) 2 ~ 20- 10 -6 and are completely negligible. This 
proves that our assumption of a non-equilibrium distribu- 
tion in the spin-degrees of freedom is justified, because 
relaxation between the hyperfine levels |5) and |6) re- 
quires a nuclear spin-flip and will therefore take place on 
a timescale set by the electron-nucleus dipolar rate which 
is much longer than the lifetime of the gas. 




B(T) 

FIG. 2. The decay rate constants as a function of magnetic 
field. Curve 1 shows G ex , whereas curve 2 gives G lai . 

In summary, we have shown that due to the large and 
negative triplet scattering length the BCS transition to a 



superfluid state occurs at experimentally accessible den- 
sities and temperatures in spin-polarized 6 Li. It should 
be pointed out however that the scattering length is ex- 
tremely sensitive to the interatomic potential, so the ex- 
act conditions required may be somewhat different than 
that given here. A better estimate of the scattering 
length can be obtained by repeating the experiment of 
Ref. || for 6 Li. Moreover, as is well-known from liquid 
3 He, the phase below the critical temperature is truly su- 
perfluid because it costs a finite amount of (free) energy 
to have gradients in the phase of the order parameter Ao. 
As a result we have a macroscopic (free) energy barrier 
for the decay of superflow and the gas can sustain persis- 
tent mass currents. By calculating the various collisional 
decay rates in the gas, we have also shown that reason- 
able lifetimes can be achieved even for densities as high as 
10 12 cm~ 3 if the bias magnetic field is larger than 10 T. 
We hope that our work will stimulate new experiments 
with this interesting quantum gas. 
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